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Magnetism in graphene is an emerging field that has received much theoretical attention. In particular, there have 
been exciting predictions for induced magnetism through proximity to a ferromagnetic insulator as well as through 
localized dopants and defects. Here, we discuss our experimental work using molecular beam epitaxy (MBE) to modify 
the surface of graphene and induce novel spin-dependent phenomena. First, we investigate the epitaxial growth the 
ferromagnetic insulator EuO on graphene and discuss possible scenarios for realizing exchange splitting and exchange 
fields by ferromagnetic insulators. Second, we investigate the properties of magnetic moments in graphene originating 
from localized -orbital defects (i.e. adsorbed hydrogen atoms). The behavior of these magnetic moments is studied 
using non-local spin transport to directly probe the spin-degree of freedom of the defect-induced states. We also report 
the presence of enhanced electron ^-factors caused by the exchange fields present in the system. Importantly, the 
exchange field is found to be highly gate dependent, with decreasing ^-factors with increasing carrier densities. 



I. INTRODUCTION 

Magnetism in graphene is an emerging field that has re- 
ceived much theoretical attention In particular, there 
have been exciting predictions for induced magnetism through 
proximity to a ferromagnetic insulator"^'^ as well as through lo- 
calized dopants and defects ^ While these systems appear to 
be very different, induced magnetic phenomena in both sys- 
tems relies on the principle of the exchange interaction. In 
the case of ferromagnetic insulators in contact with graphene, 
atomic orbital overlap at the interface is expected to induce 
a spin- splitting in the graphene layer"^"^. Alternatively, on a 
local scale, magnetism can be induced in graphene through 
dopants and defects^'^'^. In this scenario, adsorbates or lattice 
vacancies effectively remove a -orbital from the graphene 
band structure. This is known to create a localized defect state 
near the Fermi energy^' Due to the Coulomb interaction, the 
defect state is spin split leading to a spin- 1/2 populated quasi- 
localized defect state As both types of induced magnetism 
rely on the exchange interaction which requires atomic scale 
overlap of wavef unctions, materials control on this scale is ab- 
solutely necessary for realizing these new behaviors. This is 
possible through molecular beam epitaxy (MBE) which typ- 
ically has very slow deposition rates (~ A/min) allowing for 
high purity growth control in the monolayer or sub-monolayer 
range. 

Here, we discuss our work on induced spin-dependent phe- 
nomenon in graphene. First, in section II we discuss the epi- 
taxial growth of the ferromagnetic insulator EuO on graphene 
and discuss possible scenarios for realizing exchange splitting 
and exchange fields by ferromagnetic insulators. Second, in 
section III we discuss our ongoing efforts on induced mag- 
netism in graphene through localized -orbital defects. Us- 
ing non-local spin transport to probe the spin-degree of free- 
dom of the quasi-localized states, we have observed the pres- 



ence of these magnetic states. These states are coupled via 
exchange to the injected spin current which experiences an 
effective exchange field due to the presence of the magnetic 
moments. In particular, through Hanle spin precession mea- 
surements, we observe enhanced electron ^-factors caused by 
the exchange fields present in the system. Importantly, the 
exchange field observed is highly gate dependent, with de- 
creasing ^-factors with increasing carrier densities. Lastly, we 
discuss an important issue for the field of spintronics related 
to the analysis of Hanle spin precession data in the case of an 
unknown effective electron ^-factor. 

II. REALIZATION OF EPITAXIAL FERROMAGNETIC 
INSULATORS ON GRAPHENE BY REACTIVE MBE 

EuO is discussed in the literature as a possible ferromag- 
netic insulator for inducing exchange by proximity^'^'^^"^"^. 
This is understandable since EuO has excellent intrinsic prop- 
erties including large magnetization per unit volume, is highly 
resistive in the insulating state, and fits the ideal model for 
an isotropic Heisenberg ferromagnet^^. This spin splitting 
was originally predicted to be approximately 5 meV based 
on measurements by Tedrow dating back to 1986"^'^^. Very 
recently, this has been revisited by density functional calcula- 
tions that has indicated that the induced exchange splitting can 
be as large as 35 meV^. Further, there have been many fan- 
tastic theoretical suggestions for possible novel device appli- 
cations including controllable magnetoresistance"^'^^"^^, gate 
tunable manipulation of spin transport^' gate tunable ex- 
change bias^^, spin transfer torque^ Also, such exchange 
splitting is a key requirement for realizing the quantized 
anomalous Hall effect in graphene systems^^'^"^. To date, ex- 
change splitting induced by a ferromagnetic insulator in con- 
tact with graphene has not been demonstrated. 
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FIG. 1. a) RHEED image of HOPG(OOOl) substrate b) RHEED im- 
age for 5 nm EuO deposited on top of HOPG(OOOl) in the adsorption- 
limited growth regime, c) ExfoHated single layer graphene flake with 
Pd/Ti electrodes in a Hall geometry, d) Graphene resistivity, p, as a 
function of the back gate voltage. Maximum p corresponds with 
the charge neutrality point (CNP), sometimes called the Dirac point. 
Positive (negative) voltages beyond the CNP correspond to electron 
(hole) type carriers, e) Quantum Hall effect in single layer graphene 
atT=10Kand B=7 T. 



While there have been many theoretical works on global 
exchange fields in graphene due to a proximity exchange in- 
teraction caused by the presence of a ferromagnetic insulator, 
there are woefully few experimental works on this topic. Re- 
cently, the integration of EuO with sp^ carbon materials was 
demonstrated^^. Fig. 1 a) shows the reflection high energy 
electron diffraction (RHEED) image of highly oriented py- 
rolitic graphite (HOPG) substrate. HOPG consists of many 
graphene layers oriented along the (0001) axis and the sheets 
are coupled by the weak van der Waals force. Since the sub- 
strate is well oriented out of the plane, a RHEED pattern 
can be obtained despite the in-plane disorder associated with 
HOPG substrate. The surface of HOPG can be considered 
to be very similar to graphene from a growth point of view. 
Fresh surfaces of HOPG(OOl) are obtained through mechani- 
cal cleavage followed by a UHV anneal at 600 °C. The base 
pressure of the growth chamber is less than 1 x 10~^^ Torr. 
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FIG. 2. a) Gate dependent resistivity demonstrating the charging ef- 
fect at room temperature for EuO deposited on graphene FET device. 
The back gate is swept from -40 V to -1-40 V and back again at a rate 
of 0.5 V/s. b) Temperature dependence of the charging effect plot- 
ted as the relative difference in the charge neutrality point between 
up and down sweeps of the back gate voltage, c) Comparison of the 
gate dependent resistance at T=10 K for pristine and EuO deposited 
on top. d) Quantum Hall effect at T=10 K after EuO deposition. 
Black (red/grey) curve correspond to B=-i-7T (B=-7 T). 



Using a special adsorption-controlled growth regime, thin 
stoichiometric EuO films can be readily produced by reactive 
]y[gg26-29 elemental Eu flux is introduced to the HOPG 
substrate maintained at 550 °C. The elevated temperature al- 
lows for full distillation of the incident Eu atomic beam. A 
partial pressure of, P02 = 1 x 10~^ Torr, converts a por- 
tion of the Eu flux to EuO which is then thermodynamically 
favored to be deposited onto the sample surface. The re- 
maining Eu flux that is not converted, readily re-evaporates. 
In this way, the stoichiometry of the film can be maintained 
since the overpressure of Eu atoms avoids the formation of 
EU2O3 and the re-evaporation of the excess Eu avoids oxy- 
gen vacancies. Fig. 1 b) shows the RHEED pattern for EuO 
grown on HOPG(OOOl). Evident in the RHEED pattern are 
streaks with lattice constants associated with EuO(OOl) [100] 
and EuO(OOl) [110]. Thus, even though symmetry would fa- 
vor EuO(l 1 1)/HOPG(0001), growth by reactive MBE prefers 
the low energy surface plane of EuO(OOl). This is not sur- 
prising since EuO shares the rock-salt crystal structure with 
MgO, which is also known for strongly preferring out of plane 
growth along the (001) axis. 

Next, we seek to understand the impact of this high temper- 
ature growth regime on graphene field effect transistor (FET) 
devices. Graphene FET devices are fabricated in a Hall ge- 
ometry as shown in Fig. 1 c). Single layer graphene flakes are 
mechanically exfoliated onto SiO2(300 nm )/Si substrate and 
verified by Raman spectroscopy. The heavily p-doped Si sub- 
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strate acts as the back gate for tuning of the Fermi level. Pd/Ti 
(60 nm/10 nm) electrodes are fabricated using standard bi- 
layer PMMA/MMA e-beam lithography followed by e-beam 
deposition and then lift-off. Pd is a significant improvement 
over Au contacts as the surface diffusion for Pd at elevated 
temperatures is much less than Au. The Pd/Ti electrodes are 
specially capped with AI2O3. The device is characterized in 
a magneto-transport cryostat prior to EuO growth. Gate de- 
pendent resistivity measurements at T=10 K are shown in Fig. 
1 d). The charge neutrality point is Vcnp = 8 V and the av- 
erage FET mobility is 5500 cm^A^s. Pristine graphene FET 
devices exhibit minimal temperature dependence in the mini- 
mum conductivity, mobility, and CNP. Further, these devices 
exhibit the integer quantum Hall effect as shown in Fig. 1 e). 
Clear plateaus in the transverse Hall conductance in values of 
±2, ±6, and ±10 e^/h measured at T=10 K and B=7 T is evi- 
dence for 2D transport expected for single layer graphene^^'^^ 
After initial characterization, the device is loaded into the 
UHV growth chamber and ramped to 550 °C at 20 °/min. 
Once the growth temperature is reached, Eu flux is introduced 
followed quickly by P02 • EuO is allowed to grow for approx- 
imately 10 minutes (2 nm) and then the growth is stopped and 
the sample is cooled back to room temperature at 20 °/min. At 
room temperature, the sample is capped with several nanome- 
ters of e-beam MgO. The device is then removed from the 
UHV chamber, wirebonded ex-situ, and then re-loaded into 
the low-temperature, high-field magnetotransport system. In- 
terestingly, at room temperature, we observe hysteresis in the 
gate dependent resistivity (see Fig. 2 a)) which is likely due to 
a charging effect stemming from the dielectric overlay er^^'^^ . 
EuO is reported to have a static dielectric constant of 26^^. 
As shown in Fig. 2 b) the hysteresis disappears at low tem- 
perature. This is quantified by recording the difference in the 
charge neutrality point, A Vcnp, between up and down gate 
voltage sweeps. Below approximately 200 K, the difference 
is less than 1 V. Fig 2 c) compares the gate dependent resis- 
tivity at T=10 K and zero field for the pristine device and 
for 2 nm EuO deposited on top. The CNP has shifted from 
Vcnp = 8 V to Vcnp = —2 V and the average mobility 
has decreased to 3900 cm^A^s. Further, we still observe clear 
quantum Hall plateaus in the transverse conductivity, a^y, at 
the correct values expected for graphene following the equa- 
tion (Jxy = (m + ^)4e^//i, where m is the Landau level filling 
factor. The fact that the mobility decreases only slightly and 
the 2D quantum Hall effect is still observed, despite the high 
temperature deposition of EuO on top an already fabricated 
device, indicates that the growth does not drastically impact 
the graphene quality. In particular, the graphene has not func- 
tionalized into graphene oxide or some other complex struc- 
ture. 

Anomalous behavior in the Hall effect due to the overlying 
EuO material is not observed at T=10 K. Observation of such 
effects which rely on out-of-plane magnetization are likely 
hindered by the isotropic nature of the ideal Heisenberg ex- 
change which characterizes EuO. Thus, the 2 nm films used 
here, are highly dominated by shape anisotropy. While the 
results presented above for thin EuO films on graphene de- 
vices is highly promising, integration of the growth process 



with thicker dielectric films with nanofabrication procedures 
is necessary. Such advances could potentially lead to the first 
observation of the quantized anomalous hall effect in single 
and bilayer graphene^"^'^^ . 

III. MAGNETIC MOMENT FORMATION AND EXCHANGE 
FIELDS IN GRAPHENE THROUGH DOPANTS AND 
DEFECTS BY MBE 

Induced magnetic phenomena in graphene through adsor- 
bates and defects has been studied extensively theoretically 
and can be understood in terms of -orbital defects ^ These 
defects can be produced through lattice vacancies or through 
localized sp^ bonding. Experimentally, the majority of the 
progress has focused on defected graphite and graphene by 
magnetometry measurements In particular, there have 

been several reports of ferromagnetic ordering in graphene 
and graphite^^""^^. However, theoretically, the ferromagnetic 
ordering in realistic experiments is not well understood^ 
Magnetometry measurements have the advantage of directly 
measuring the total magnetic moment, but are subject to ar- 
tifacts. The recent work of Nair, et al.,^ demonstrated that 
many HOPG substrates, from which graphene flakes are de- 
rived from, are littered with magnetic impurity clusters which 
could possibly explain the reports of ferromagnetic behavior 
in magnetometry measurements. The authors went on to con- 
clusively demonstrate spin, S=l/2, paramagnetism induced in 
graphene laminates through proton irradiation and fluorina- 
tion. Alternatively, there have been several transport measure- 
ments on doped and defected graphene. This has the advan- 
tage of locally probing the graphene flake, but does not di- 
rectly probe the magnetic state in the way that magnetometry 
does. In 2011, Candini et al.,"^^ observed hysteretic magne- 
toresistance in narrow graphene ribbons and Hong et al.,^^ 
observed negative colossal magnetoresistance and localiza- 
tion in dilute fluorinated graphene. Notably, Chen et al.,"^^ 
reported the observation of gate tunable Kondo effect in de- 
fected graphene through measurements of the temperature de- 
pendent resistance which demonstrated a possible logarith- 
mic dependence that could be suggestive of the Kondo effect. 
This claim was disputed by Jobst and Weber"^^, and Kondo 
in graphene remains a controversial topic. Very recently, we 
have demonstrated magnetic moment formation in doped and 
defected exfoliated graphene flakes by combining the advan- 
tages of these two approaches^. More specifically, we seek to 
probe the spin degree of freedom locally to the graphene in 
order to minimize the potential for artifacts and minimize the 
uncertainty in interpretation of the experimental data. 

The proposed experiment is schematically displayed in Fig. 
3 a). Non-local spin valves are employed to directly probe the 
spin degree of freedom locally on the graphene flake. The 
devices can then be loaded into a UHV deposition chamber 
(base pressure 1 x 10~^^ Torr). The chamber is specially 
equipped with a sample manipulation stage that is attached 
to a cryostat. Electrical probes can be connected to the de- 
vice inside the chamber without exposing the sample to air or 
needing to vent to manipulate the electrical probes. There is 
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FIG. 3. a) Schematic of the proposed experiment in which atomic 
hydrogen is carefully introduced in-situ to non-local graphene spin 
valves in an MBE chamber with magnetotransport capability b) Gate 
dependent resistivity of a pristine non-local graphene spin valve be- 
fore the introduction of atomic hydrogen, c) Non-local resistance, 
i^A^L, as a function of in-plane applied field, Bapp,y, along the axis 
of the magnetic electrodes for pristine graphene. The two state par- 
allel and antiparallel behavior demonstrates spin injection and trans- 
port in the graphene channel. Data corresponds with the channel in 
b) maintained at Vg = —15 V. d) Hanle spin precession measure- 
ments in pristine graphene at Vg = —15 V. d) Gate dependence of 
the spin lifetime determined from fits to the spin precession data for 
pristine graphene with tunneling contacts using the Hanle equation. 
All data is taken at T=12 K. 



also a home-made Fe core magnet inside the chamber. Thus, 
magnetotransport measurements on non-local graphene sam- 
ples can be performed at cryogenic temperatures immediately 
before and after -orbital defects in the form of atomic hy- 
drogen is introduced to the graphene surface. The entire ex- 
periment can be performed without needing to warm up or 
expose the sample to air. Atomic hydrogen is produced from 
H2 using an Omicron source. High voltage deflectors lo- 
cated in front of the atomic hydrogen source steer charged 
ions away from the sample. Fabrication and characterization 
of graphene non-local spin valves with tunneling contacts is 
reported in detail elsewhere^'"^^"^^. Tunneling spin valves, or 
high resistance contacts, are crucial for this experiment as 
they alleviate the conductivity mismatch problem and increase 
the observed spin lifetime in graphene by nearly one order of 
magnitude"^^'"^^'^^. 

Fig. 3 b) displays the gate dependent resistivity for a 
graphene spin valve measured in the four terminal geometry 



using standard lock-in techniques. The non-local resistance is 
measured in the configuration displayed in Fig. 3 a). Elec- 
trodes B and C are ferromagnetic Co with MgO tunnel barri- 
ers, while A and D are spin insensitive Au contacts. Current is 
passed from electrode B to A, with electrode A held grounded. 
Spins are injected into the graphene under electrode B. A spa- 
tial gradient in the chemical potential causes spins to diffuse 
to the right along the channel towards electrode C where they 
can be detected as a voltage relative to electrode D. A non- 
local resistance, Rnl can be defined as the non-local voltage, 
Vatl, measured between electrodes C and D, divided by the 
applied current (RNL=yNL/^)- Sweeps of the in-plane mag- 
netic field along the electrode length generates characteristic 
two state Rnl curves as shown in Fig. 3 c) measured at T=12 
Kand Vg = -15 V. 

Instead of applying an in-plane field, we can apply an out- 
of-plane magnetic field, Bapp^z^ in which case the injected 
spins (oriented in-plane collinear to the magnetization axis of 
the electrode wires) precess around the applied field. The dif- 
fusive transport causes different arrival times at electrode C. 
When measured over time scales longer than the spin relax- 
ation and transit time, this appears as dephasing of the spin 
population with increasing applied field. A typical spin pre- 
cession measurement taken at T=12 K and Vg=-15 V is shown 
in Fig. 3 d) for parallel (blue) and antiparallel (red) relative 
magnetization alignments. For high resistance contacts, the 
spin precession measurements on pristine graphene are well 
explained by the Hanle equation^'"^^. 



Rnl = Snl 



^-L'^/ADt 

V^TrDt 



cos(^^ Bapp,zt)e-'^^^dt (1) 



where is of spin lifetime, D is the diffusion coefficient, 
Sn l is the Hanle amplitude, A = ^/Dts is the spin diffusion 
length and g is the electron ^-factor. For pristine graphene, 
g is assumed to be 2. Thus, by fixing the electron ^-factor, 
which is well known for pristine graphene^ \ the Hanle data 
can be uniquely fit using equation 1 and values for the spin- 
lifetime and the spin diffusion length can be obtained. In this 
case, the diffusion coefficient is renamed D=Ds to signify 
that it is determined from the spin measurement. Analysis of 
the data presented in Fig. 3 d) yields rs=479 ps, D=Ds=0.023 
m^/s and A=3.3 fim for the channel length, L=5.25 /im. The 
spin lifetime for the pristine graphene device is plotted for sev- 
eral gate voltages in Fig. 3 e) and best fits yield spin lifetimes 
in the range of 300 to 600 ps. Here we distinguish Ds from 
Dc, the diffusion coefficient associated with the charge trans- 
port. Dc can be experimentally determined using the Ein- 
stein relation, Dc=cf /e^v, where cr is the measured channel 
conductivity, e is the electron charge, and v is the density of 
states (DOS) for single layer graphene. For pristine graphene, 
Dc and Ds values are expected to agree^'^^ 

Next, we examine the effect of atomic hydrogen on charge 
and spin transport. Atomic hydrogen doping and all measure- 
ments occur at T=12 K. Fig. 4 a) displays the gate dependent 
resistivity for several atomic hydrogen exposures. The blue 
curve is the pristine device (0 s) and corresponds with the data 
in Fig. 3 b). The resistivity dramatically increases with 2 s 
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and then again with 8 s hydrogen exposure. The charge neu- 
traUty point shifts minimally from Vcatp^O V to VcArp=-l 
V and the minimum conductivity changes dramatically. This 
is characteristic of the formation of -orbital defects as op- 
posed to charge impurity (CI) scattering^^'^^ . From the reis- 
tivity, we can obtain an order of magnitude estimate for the 
number of -orbital defects of 0.1%. Accordingly, the mo- 
bility is dramatically decreased when the band structure is so 
affected locally as shown in Fig. 4 b). The decrease in the 
mobility accompanies a decrease in the diffusion coefficient, 
which is expected to decrease the overall non-local resistance 
signal. Thus, as shown in Fig. 4 c) the non-local resistance, 
Rnl^ decreases with increasing hydrogen exposure for Vg=- 
16 V. The key feature in \hQ Rnl data is the decrease in sig- 
nal around zero field, called the dip in Rnl • With increasing 
atomic hydrogen exposure, the size of the dip increases rela- 
tive to the size of the non-local signal at high fields. 

This dip is due to the creation of localized magnetic mo- 
ments that can interact with the injected spin current via 
exchange coupling^. A model Hamiltonian including spin- 
moment interactions for a single spin interacting with an en- 
semble of induced moments can be written as. 



H = r]MAexSe • (Sm) + gj^sSe ' Ba 



(2) 



where rfu is the filling density of magnetic moments, Aex is 
the exchange coupling between the conduction electron spin 
(5'e) and the spin of the induced moment {Sm)- The effec- 
tive field generated by the exchange interaction is given by 
Bex = "^mA^A^^ In the local frame of the electrons, the ex- 
change field can be described by a time-dependent, randomly 
fluctuating magnetic field, Bex{t), that is the sum of an aver- 
age mean field term and a randomly fluctuating component. 



Bex{t) = Bex + ^Bex{t). 



(3) 



These two components of the exchange field impact the ob- 
served data in different ways. The exact line-shape of Rnl 
vs. Bapp^y and the dip are explained by the randomly fluctu- 
ating effective exchange field. This fluctuating effective field 
is characterized by net-zero average value and non-zero rms 
fluctuations. This is known to generate extra spin relaxation 
that depends in a characteristic fashion (Lorentzian) upon the 
applied field {Bapp^yf'^^ . On the other hand, the mean field 
term creates rapid Larmor precession in Hanle measurements, 
and will be discussed in more detail later in the article. In Fig. 
4 d) we show the gate dependence of the relative dip size for 2 
s and 8 s hydrogen exposure. This is defined quantitatively 
as RNL{Bapp,y = 0)/RNL{Bapp,y ^ oo) and 100% indi- 
cates complete scattering (relaxation) of the spin population 
and zero spin signal at zero field. The effect on the spin scat- 
tering at zero field is symmetric for electrons and holes and 
the dip size decreases with increasing carrier concentration. 

We now turn our attention to Hanle spin precession mea- 
surements for hydrogen doped graphene. The Hanle curves 
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FIG. 4. a) Gate dependent resistivity curves for pristine graphene 
(blue) and graphene exposed to 2 s (black) and 8 s (red) of atomic 
hydrogen. The pristine curve (blue) corresponds to the data pre- 
sented in Fig. lb), b) The electron (black) and hole (red) mo- 
bilities for the graphene spin valve as a function of hydrogen ex- 
posure, c) Non-local resistance, Rnl = —16 V for 2 s 
(left axis, black) and 8 s (right axis, red) of atomic hydrogen ex- 
posure. Notably, the dip at zero field increases with increasing hy- 
drogen exposure, d) Gate dependence of the relative dip size com- 
pared to the high field Rnl values. This is defined quantitatively 
as RNL{Bapp,y = 0) / Rn L{Bapp,y oo) . 100% corrcsponds with 
full spin scattering and zero spin signal at zero field. 



at Vg=-16 V, -1 V, and 14 V are displayed in Fig. 6 a), b), 
and c), respectively. For pristine graphene (see Fig. 3 d)), 
a precession of 180° occurs between 30 mT and 40 mT for 
all gate voltages. On the other hand, for hydrogen-doped 
graphene, the Hanle curves are dramatically sharpened with 
a full 180° rotation occurring within 5 mT at all gates. At first 
glance, the sharpened Hanle is suggestive that the exchange 
field is playing a role in the observed spin precession. We 
now run into a key issue with Hanle analysis. Standard Hanle 
fitting employs equation 1 and assumes g=2 while allowing 
the amplitude, lifetime, and diffusion coefficient to be fitting 
parameters. However, with the possibility that g ^ 2, then 
we must replace g with an effective ^-factor, This is prob- 
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lematic for Hanle fitting since non-unique fits can be obtained 
with Snl, D, Ts, and ^* all free parameters. This can be 
seen explicitly if we examine the transformation, g"^ ^ cg"^, 
T2 ^ T2/ c, Ds c Ds, S ^ cS, where c is a constant. For 
a given parameter set (^* , T2 , Ds, S), the transformed param- 
eter set {eg* ^T2/c^cDs,cS) will generate the same Hanle 
curve. Therefore, if we force a Hanle fit that assumes g=2, a 
transformation that effectively decreases ^* (i.e. c<l) has the 
effect of artificially enhancing the resulting value and arti- 
ficially decreasing the Hanle diffusion coefficient, Ds- This 
effect is evident in Fig. 5 a) in which the diffusion coeffi- 
cient for 8 s hydrogen exposure is plotted for both Dc and 
Ds- Here, Ds is obtained from fitting the 8 s hydrogen-doped 
Hanle data and fixing g = 2, while Dc is obtained from the 
measured graphene conductivity using the Einstein relation. 



Dc = 



(4) 



Due to charge inhomogeneities in the Si02/Si substrate, 
we must use the broadened DOS, y{E). For Gaussian 
broadening^^, the DOS is 



v{E) 



2b 



27r 



(5) 



where gy is the valley degeneracy, gs is the electron spin de- 
generacy, h is Planck's constant, Vi?=l x 10^ m/s is the Fermi 
velocity, and 6=100 meV is the Gaussian broadening parame- 
ter. Accordingly, the calculated Dc shows reasonable agree- 
ment with the conductivity (see Fig. 4 a). Unlike the case 
of pristine graphene^'^^, Ds and Dc are drastically different 
as summarized in Fig. 5 b). Such a discrepancy is a key in- 
dicator that ^* 7^ 2. This an important issue for spintronics 
as standard Hanle fitting procedures assume the material g- 
factor to be fixed. This is directly relevant to recent exper- 
iments concerning -orbital defected graphene^^'^^, but has 
also been discussed for epitaxial graphene on SiC^^'^^. For 
SiC, the graphene channel may be strongly coupled to local- 
ized states in the underlying substrate^^. In that system, there 
exists a large discrepancy between Dc and Ds, and enhanced 
effective ^-factors are observed. 

In accordance with the above discussion and following the 
procedure outlined by Maassen et al , we obtain fits to the 
Hanle spin precession measurements using equation 1 with 
free parameters Snl^'^s^ and Dc is fixed to be the exper- 
imentally determined value. The fits (red curves) are shown 
in Fig. 6 a), b), and c) and describe the observed Hanle spin 
precession data extremely well. In Fig. 6 d) we plot the re- 
sulting ^* values as a function of gate voltage for 2 s (blue 
diamonds) and 8 s (red triangles) hydrogen-doped graphene. 
The ^* gate dependence is strikingly similar to the relative 
dip size (see Fig. 4 d)) indicating that the two phenomena are 
indeed related. This agrees well with the theoretical picture 
proposed above. With increasing carrier concentration, the ef- 
fective exchange field, ^e^, = ^^^^^"^'^^^ = B ex exit), 
decreases. This impacts both the non-local resistance in plane 
field measurement and the Hanle spin precession measure- 
ments. While the nature of the gate dependence remains an 
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FIG. 5. a) The experimentally determined diffusion coefficient from 
the Einstein relation (Dc) (black curve) and the diffusion coefficient 
obtained from Hanle fitting with g=2 (Ds) (black squares) plotted as 
a function of back gate voltage for 8 s hydrogen exposure. For com- 
parison, the measured conductivity (blue curve) is shown on the right 
axis, b) Relative difference between Ds and Dc for 8 s hydrogen- 
doped (dark red squares) graphene. 
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FIG. 6. a) - c) Hanle spin precession curves for graphene spin valves 
exposed to 8 s atomic hydrogen at several gate voltages. Black open 
circles represent the measured Rnl signal vs. Bapp,z. Red line is 
the fit to the Hanle equation using fixed Dc (determined from the 
experimentally measured resistance). The spin lifetime and electron 
g-f actor allowed to be free parameters, d) ^* values determined from 
Hanle fitting for 2 s (blue diamonds) and 8 s (red triangles) atomic 
hydrogen exposure. Dashed line indicates g=2. e) Spin lifetime ob- 
tained from Hanle fitting for pristine graphene (black squares), 2 s 
hydrogen doping (blue diamonds), 8 s hydrogen doping (red trian- 
gles). The pristine data is the same as shown in Fig. 3 e) with g=2. 
The spin lifetimes after hydrogen doping for 2 s and 8 s corresponds 
with the g* values shown in d). 



open question, the above equation suggests several factors that 
could play a role. It is possible that as the Fermi level is tuned 
away from charge neutrality, that either the -orbital defect 
states become double occupied for positive gates or unoccu- 
pied with negative gates. In either case, this would eliminate 
the local spin moments. This scenario has been supported very 
recently by magnetometry experiments'^^. Alternatively, the 
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exchange coupling may be strongly affected by hybridization 
or screening, which may further decrease the signal at high 
carrier concentrations. Here screening may play a significant 
role by reducing the exchange overlap between a carrier spin 
traversing a region in the channel and the induced moments at 
further distances. 

The spin lifetime is also obtained through this fitting proce- 
dure. Fig. 6 e) displays Ts for 2 s (blue diamonds) and 8 s (red 
triangles) hydrogen exposure. For comparison, the pristine 
graphene lifetime is displayed (black squares). Notably, there 
is a slight trend upwards in spin lifetime of roughly a factor 
of 2 with up to 8 s hydrogen exposure. The average spin life- 
times are 377 ± 69 ps, 548 ± 96 ps, and 777 ± 130 ps, over 
this back gate range for pristine, 2 s, and 8 s, respectively. 
This may indicate the importance of the D'yakonov-Perel^"^'^^ 
spin relaxation mechanism in single layer graphene. Unfor- 
tunately, this cannot be considered conclusive evidence as the 
exchange field complicates the analysis. Recent experiments 
have investigated spin relaxation mechanisms in high mobility 
graphene^^'^^ and charge impurity doped graphene^^'^^. How- 
ever, despite the ongoing efforts, this issue remains an open 
question in the field of graphene spintronics. 

Lastly, we note that this analysis follows a completely dif- 
ferent procedure than we employed in McCreary et ah ^ on 
the same Hanle data. In the previous analysis, the Hanle fit- 
ting utilized parameters determined from a unique fitting pro- 
cedure based on analysis of the dip observed in the in-plane 
Rnl scans as well as involving data taken prior to hydrogen 
doping. The current analysis is much simpler and does not 
rely on in-plane Rnl scans or data taken prior to hydrogen 
doping. Importantly, both methods of Hanle analysis yield the 
same result that ^* is much greater than 2 and has a maximum 
value near the Dirac point. Because the current Hanle anal- 
ysis does not involve the in-plane Rnl scan, the observation 
of enhanced ^* provides direct evidence for the formation of 
magnetic moments that is completely independent of the zero 
field dip in the in-plane Rnl scan. 

IV. CONCLUSION 

Induced spin-based phenomena in graphene is an emerg- 
ing field with great potential for realizing new physics and 
device structures. Here we have investigated induced mag- 
netism in graphene and the exchange interaction, which relies 
on wavefunction overlap. For this, we have employed MBE 
in order to realize highly controlled interfaces and materials. 
First, we have examined the integration of EuO onto graphene 
FET devices in a high temperature adsorb tion-limited regime. 
It was found that at low temperatures the device is not dramat- 
ically impacted and still behaves nearly as well as the pristine 
graphene device. This could lead to the future realization of 
the exchange proximity effect over large regions of a graphene 
layer. Further, we have examined induced magnetism on a lo- 
cal scale through the systematic introduction of -orbital de- 
fects in the form of atomic hydrogen. By employing non-local 
spin valves in a special in-situ magnetotransport cryostat and 
MBE growth chamber, we have been able to systematically 



investigate the interaction of pure spin currents with the in- 
duced moments in a highly controlled environment, without 
ever needing to warm up the sample or expose to air. Here, 
the carrier spins interact with the induced moment spin state 
through the exchange coupling. Effects of this coupling are 
observed in both in-plane non-local resistance scans and in the 
Hanle spin precession data. In particular, we have investigated 
the mean-field component of the exchange field which leads 
to enhanced effective electron ^f-factors and fast spin preces- 
sion. The enhanced ^-factors are gate dependent exhibiting a 
symmetry around the charge neutrality point with decreasing 
precession frequency with increasing carrier concentrations. 
This could possibly be explained by changes in the occupancy 
of the defect states as the Fermi level is changed or through a 
screening effect which decreases the coupling between car- 
rier spins and the induced moments. Lastly, here we have 
proposed a simpler Hanle analysis method compared to a pre- 
vious treatment of exchange coupling to induced moments in 
graphene spin valves by McCreary et al ^ . Notably, this sim- 
pler procedure produces the same physical picture and also 
yields enhanced electron ^-factors in the range of 5 to 15, well 
away from the pristine value of 2. These results demonstrate 
the importance of MBE for investigating spin-based phenom- 
ena on the atomic scale. 
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